Higher-dimensional chaotic stadium billiards with cylindrical shape 
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We describe conditions under which higher-dimensional bilhard models in bounded, convex regions are 
fully chaotic. Our models generalize the Bunimovich stadium to dimensions above two. An example is a 
three-dimensional stadium bounded by a cylinder and several planes, the combination of which give rise 
to a defocusing mechanism. We provide strong numerical evidence that this and other such billiards are 
fully hyperbolic — all but two of their Lyapunov exponents are non-zero. Applications to tubular networks 
and other cavities, as well as to models of interacting particles, are discussed. 
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Billiard models, in which a point particle moves freely 
between elastic collisions with a fixed boundary, are a fer- 
tile source of ideas in physics [IJ and mathematics lEl 
alike. They provide a basis for some of the fundamental 
concepts of statistical mechanics, and are at the same time 
open to mathematically rigorous analysis [3 1. In particular, 
they are some of the best-motivated models which exhibit 
strong chaotic dynamics. 

There are two distinct categories of chaotic billiards. The 
first is the class of dispersing billiards, the prototypical ex- 
ample of which is the hard-sphere gas: the dynamics of N 
hard spheres in a three-dimensional box with elastic colli- 
sions is equivalent to a point particle in a 3A/^-dimensional 
space moving uniformly outside a collection of spheri- 
cal cylinders, with specular reflections at the boundary 
Two-dimensional examples include the Sinai billiard 
lO and the periodic Lorentz gas [6|, in which a particle 
bounces off a disk on the 2-torus or a periodic configu- 
ration of them on the plane. The latter has fast decay of 
correlations and thus converges to Brownian motion, i.e., 
is diffusive in a strong sense [7 J. The mechanism giving 
rise to chaos in such billiards is that of dispersion, where 
nearby trajectories separate at each collision with a convex 
surface; this leads to an overall exponential divergence, or 
equivalently to positive Lyapunov exponents, and the sys- 
tem is then said to be completely chaotic or hyperbolic [8 |. 

The second category is made up of defocusing billiards, 
the most well-known example of which is the Bunimovich 
stadium [9|. Here, chaos is due to a mechanism different 
from dispersion — that of defocusing. Unlike in the Sinai 
billiard, the boundary of the stadium curves inwards with 
respect to the particle. Nearby trajectories initially focus 
after colliding with this boundary; however, the distance to 
the next collision is typically longer than the distance to the 
focal point, so that they eventually defocus even more. This 
again leads to an overall exponential expansion in phase 
space and hence complete chaoticity. 

Defocusing billiards have attracted much attention in the 
physics community, particularly in connection with quan- 



tum chaos fWl. The stadium billiard has served as a 
paradigm for determining, both theoretically and exper- 
imentally, the statistics of the eigenvalues and eigenvec- 
tors of classically chaotic systems [TTl [T2l[T3l[T4l[T5l[T6ll , 
and for investigating quantum localization |[T71[T8l. Stadia 
have also been studied in acoustic experiments in closed 
chaotic cavities fil9l and optical microcavity laser exper- 
iments EOl EB El El El; see in particular Ref. [25] 
for three-dimensional cavity experiments. Stadium- shaped 
microstructures have been used in quantum conductance 
experiments ll26l l27l [28l , and ray dynamics were studied 
in Ref. fM . 

The extent to which the defocusing mechanism works 
in dimensions beyond two has, however, long remained 
unclear. Examples of three- and four-dimensional chaotic 
billiards with flat and spherical components were proposed 
and numerically investigated in [30|; see also [31 1. It was 
subsequently proved that defocusing can produce chaos in 
higher-dimensional billiards [32], but the examples given 
were non-convex. A particular three-dimensional convex 
stadium billiard, with two perpendicular half-cylindrical 
components was numerically investigated and found to be 
chaotic in Ref. [i33J , and this was recently proved rigor- 
ously in Ref. [|34l . Its construction is, however, compli- 
cated and rather non-generic. 

In this Letter, we establish the existence of a large class 
of higher-dimensional, convex billiards which are defocus- 
ing and fully chaotic. They are based on cylindrically- 
shaped structures, which, as we argue below, give rise to 
chaos in conjunction with planar elements. Our results 
demonstrate that hyperbolicity in defocusing billiards is 
easier to obtain than was previously believed. 

The construction of our class of models is simplicity it- 
self: in dimension three, they are formed by cutting a cylin- 
der with one or more flat planes to form a convex region, 
inside which the billiard dynamics takes place. Simple ex- 
amples, consisting of a three-dimensional circular cylinder 
cut by two and three planes, are shown in Fig. [T] In both 
cases, one of the planes is perpendicular to the cylinder 
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FIG. 1 : (Color online) Three-dimensional stadia consisting of a 
circular cylinder cut by (a) two planes, one perpendicular to the 
cylinder axis and the other at angle 7r/4, and (b) three planes, 
one perpendicular and two intersecting at right angles and cutting 
the cylinder at angle 7r/4. In each case, a typical trajectory is 
depicted, with varying shades as time progresses. 

axis, which just serves to confine the motion without in- 
troducing any new dynamical features. The other planes, 
however, are angled away from perpendicular. This is al- 
ready sufficient to render the system completely chaotic. 

One might naively expect a cylindrical- shaped billiard 
to produce only integrable behavior. This would indeed be 
the case if the planes were all perpendicular to the cylin- 
der axis, leaving invariant the angular momentum along the 
cylinder axis. However, when one of the planes is oblique 
with respect to the cylinder axis, integrability is lost and 
defocusing can take place. 

In general, we call a cylindrical stadium billiard a 
bounded, convex region made by cutting a cylinder with 
flat planes, such that at least part of the boundary of the 
region is curved, and such that the symmetries of the sys- 
tem are broken. This construction easily extends to higher 
dimensions, as discussed below. 

We conjecture that generically the classical billiard dy- 
namics within such a region is completely chaotic. That is, 
all of its Lyapunov exponents, which measure the separa- 
tion rate of nearby trajectories in phase space, are nonzero, 
with the exception of the two associated with energy con- 
servation and the corresponding time translation. This 
claim is substantiated by the results of numerical computa- 
tions of these exponents, such as those shown in Fig.|2] 

The mechanism leading to hyperbolicity in these mod- 
els is a form of defocusing. The oblique planes allow 
the expansion of wavefronts emanating from the cylindri- 
cal surface before they recollide with that surface. This 
can be visualized by "unfolding" the billiard, i.e., reflect- 
ing it in its planes. For example, unfolding the billiards 
of Fig. [T] gives equivalent cylindrical square- and cross- 
shaped billiards, respectively, as shown in Fig. [3] The 
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FIG. 2: (Color online) Spectrum of Lyapunov exponents of the 
billiard shown in Fig.[TJb). The parameter p measures the ratio 
between the radius of the cylinder and half its height. There are 
four non-zero exponents, arranged in two positive-negative pairs. 
The inset shows a log-log plot of the positive exponents. As p 
increases, the exponents decrease with the frequency of collisions 
with the oblique planes. 

latter can be viewed as a building block to generate ex- 
tended arrays of shuch a cylindrical structure, by reflecting 
it in its flat planes. Within a cylinder, the dynamics is of 
the integrable focusing-defocusing type, such as found in 
a two-dimensional circular billiard. When trajectories are 
reflected in the angled plane, however — or, equivalently, 
continue on to the next collision when unfolded — nearby 
trajectories can have time to defocus before colliding again 
with the cylinder surface. 

We emphasize that our constructions are truly higher- 
dimensional. Apart from Ref. (301, which considered a 
sphere cut by planes, previous constructions have mainly 
consisted of twisted Cartesian products [ 35 1. 

We expect that a rigorous proof of the hyperbolicity, and 
indeed ergodicity, of our models is within reach by exploit- 
ing this idea, using the techniques of Refs. |[34l|35j|. This 
is however not immediate as our models do not satisfy the 
separation conditions required for the proofs in Ref. [^1: 
there is no spatial separation between parts of the cylin- 
drical surface, nor is there separation in time, since a tra- 
jectory can bounce near the intersection with a plane and 
return arbitrarily quickly to collide again with the cylinder 
surface [40|. 

Examples of the general class of billiards described 
above arise naturally in models of interacting particles, 
such as previously discussed in |[37ll . The addition of flat 
surfaces, however makes them somewhat more general. 

Consider, for example, two point particles in two di- 
mensions, joined by a massless string of length a and 
trapped between parallel walls separated by a distance d. 
The particles interact only when the string becomes taut, 
at which point they exchange momentum parallel to the 
string. Denoting by x/ = the particle positions, the 

dynamics is confined to the region given by the inequali- 
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FIG. 3: (Color online) Cylindrical (a) square- and (b) cross- 
shaped billiards obtained by unfolding the billiards shown in 
Fig. [T] with representative chaotic trajectories. Other similar 
structures, including networks with joins, can be obtained by cut- 
ting the cylinders with planes at different angles. 

ties ||xi — X2IP < and \yi\ < d/2 (i = 1,2). The transla- 
tional symmetry of the system implies conservation of the 
momentum of the center of mass parallel to the channel 
walls. Considering only the relative motion of the two par- 
ticles between the planes, the phase-space dimensionality 
reduces to six, and the dynamics is equivalent to billiard 
dynamics in a three-dimensional convex region, given by 
x^+y^ < {a/d)^ and < 1, where the coordinates 

are given by x= {xi —X2)/d, y = (yi —y2)/d and 
z = {yi +y2)/d. This billiard is equivalent to that shown in 
Fig.[TJb), after reflection in the bottom plane. 

The generalization to higher-dimensional billiards is 
straightforward. In four dimensions, for example, there 
are three families of hypercylinders. We could thus have 



billiards made either of a spherical cylinder (a three- 
dimensional spherical base with a single axis), or of a cylin- 
drical prism (a two-dimensional circular base and two per- 
pendicular axes), both cut by flat (hyper)planes breaking 
the symmetries of the cylinders. The third family corre- 
sponds to billiards made of the Cartesian product of two 
disks, which are in fact equivalent to the model of Ref . f3T| 
with three particles in two space dimensions, after the di- 
mensions associated to three of its invariants (the angular 
momentum and the two components of the momentum) 
have been removed. In each of these cases we conjecture 
that the dynamics will generically be chaotic (with three 
pairs of positive-negative Lyapunov exponents). 

An example of such a higher-dimensional billiard is 
given by confining the interacting particles used above to a 
square box; see Fig.|4] This gives rise to a four-dimensional 
cylindrical prism billiard, with each end cut by four planes, 
each at angle n/A to an axis of the cylinder. Numerical 
calculations show that the system is again fully chaotic. 

Similar results are obtained for two particles in a three- 
dimensional box, which is equivalent to a six-dimensional 
hypercylinder with a three-dimensional spherical base and 
three perpendicular axes. Chains of multiple particles con- 
fined to boxes, such as those investigated in the context of 
heat conduction in Ref. |[38]| . also yield higher-dimensional 
chaotic stadia. 




FIG. 4: (Color online) Dynamics of two particles in a square box. 
Two nearby trajectories are shown to diverge quickly after collid- 
ing with the walls. The solid lines show the particle trajecto- 
ries, and the dashed lines the centers of mass. The corresponding 
billiard is a four-dimensional hypercylinder of cylindrical prism 
type, cut by pairs of intersecting hyperplanes. 

A main conclusion of our work is that the integrability of 
billiards in two dimensions can be easily destroyed in the 
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presence of an additional spatial dimension, provided that 
symmetry is broken. Using the basic idea that breaking 
cylindrical symmetries is enough to induce chaos, we have 
shown that the class of completely chaotic, convex, three- 
dimensional billiard models, or stadia, is much larger than 
previously believed. 

Although the cylinders used here have circular bases, we 
expect that cylinders with other convex bases, such as two- 
dimensional stadia, are also chaotic. Furthermore, prelim- 
inary results suggest that cylinders with other sufficiently 
smooth, convex bases, such as ellipses, can also be com- 
pletely chaotic. Extra dimensions are thus also able to de- 
stroy elliptical islands present in the base system. This will 
be the subject of future work. 

We have also shown that our billiards can form build- 
ing blocks for spatially-extended periodic and nonperiodic 
structures; these are related to the track billiards introduced 
in Ref. [39], but with piecewise straight segments, angled 
with respect to one another, and providing a mechanism 
for particles to turn around. Our results imply that classical 
dynamics within such structures is chaotic. This can be ex- 
pected, for instance, in a series of straight tubes connected 
by joins. Such extended structures display diffusive behav- 
ior, as will be reported elsewhere. The existence of such 
chaotic structures has further interesting implications for 
physical systems and many potential applications, whether 
in nano-devices, fluid dynamics, acoustic devices or optical 
fibers, where experiments would be possible. 
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